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ABSTRACT 

A  numerical  method  for  solving  a  problem  In  unsteady  slag  flow  in  the  /  t>f'c 

(  t'*'1 

hearth  of  a  blast  furnace  is  presented.  This  problem  is  reduced  to  a  free  \ 
boundary  problem  for  an  elliptic  system.  The  potential  problem  for  a  given 
free  boundary  is  approximated  by  the  penalty  method.  The  derivatives  of  the 
potential  function  on  the  free  boundary  is  approximated  by  the  integration  of 


the  penalty  term,  and  then  the  subsequent  shape  of  the  free  boundary  is 
obtained  by  solving  the  differential  equation  for  the  motion  of  the  free 
boundary.  The  finite  difference  method  is  used  to  solve  the  penalized 
problem.  A  numerical  example  is  given.  <^-' 
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NUMERICAL  SOLUTION  OF  FREE  BOUNDARY  PROBLEM 
FOR  UNSTEADY  SLAG  FLOW  IN  THE  HEARTH 


•  •* 

Makoto  Natori  and  Hideo  Ka war a da 

1.  INTRODUCTION 

We  present  a  numerical  method  for  solving  a  problem  in  unsteady  flow  of 
molten  slag  in  the  hearth  region  of  iron  producing  blast  furnaces  during  the 
tapping  operation  [ 1 ] .  This  problem  is  reduced  to  a  free  boundary  problems 
for  an  elliptic  system.  This  type  of  problem  is  similar  to  the  porous  flow  of 
underground  water  in  which  the  water  surface  is  a  free  boundary.  The 
numerical  calculations  of  this  type  of  problem  were  done  by  various 
researchers  [2  -  5] .  The  three-dimensional  problem  of  the  slag  flow  in  the 
hearth  was  solved  by  using  the  finite  element  method  by  Ichihara  and  Fukutake 
[6].  They  concluded  that  their  computation  scheme  is  not  efficient  in 
practical  use. 

The  object  of  this  paper  is  to  resolve  this  computational  instability  by 
using  the  penalty  method  developed  by  Kawarada  and  Natori  [7  -  10]. 

2.  FORMULATION 

We  consider  two-dimensional  slag  flow  in  the  hearth  which  is  bounded  by 
impermeable  boundaries  y  *  0,  x  *  0  and  x  **  a.  One  of  vertical  boundaries* 
x  *  0,  has  a  tapping  hole  near  the  bottom.  As  shown  in  Figure  1,  y  *  g(x,t) 
denotes  the  free  surface  of  the  slag  region  0^: 
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n  -  { (x,y)  |  0  <  X  <  a,  0<y<  g(x,t)} 
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Free  boundary 
y  =  g(x,t) 


Tapping  hole 


a  x 


Figure  1 


A  velocity  potential  can  be  defined  by 

♦  *  ^  +  y 

Y  1 

where  p  is  the  fluid  pressure  and  Y  is  the  specific  weight  of  the  fluid 
If  it  is  assumed  that  Darcy's  law  holds,  the  potential  is  given  by 

(1)  A  *  -  0  in  JJ 

g 

(2)  $  *  y  on  y  *  g(x,t) 

(3)  $  »  0  on  y  «  0 

y 

(4)  -  0  on  x  ■  0  and  x  -  a,  except  on  the  tapping  hole 


$x  *  k  (  >0  )  on  the  tapping  hole 


where  k  is  a  constant. 


V 


defined  later 


Figure  2 


e 

By  the  use  of  x  ,  equation  (1)  is  approximated  by 

(8)  A*e  -  -  x€($e  '  y)  ■  0  in  0  , 

where  £  is  a  positive  constant.  We  add  a  new  boundary  condition: 

e 

(9)  +  -  y  on  y  *  b  , 

to  the  boundary  conditions  (3)  -  (5). 

e 

In  fact,  if  we  let  e  be  suffciently  small  then  we  know  that  + 

e 

approximates  ♦  in  8  c  and  4  is  nearly  equal  to  y  in  8-8  e  [12]. 

9  9 

Therefore  the  boundary  condition  (2)  is  approximately  satisfied. 

If  we  use  the  method  of  integrated  penalty,  equation  (6)  is  approximated 
by 

(10)  g*  -  -d  -  ~  Jq  X6(4E  -  y>dy>  • 

This  equation  is  obtained  as  follows.  We  put 


e 

P 


* 


e.  i  1  e..e  . 
"  p  (x,y)  -  “  X  (♦  -  y) 


<1  m  qC(x,y)  -  f*  p6(x,n)dn  . 

By  an  application  of  Theorems  1.1  and  1.2  in  [11],  we  have 


»'  *  £  in  • 


*-  Si  1/1  ♦  4 11  -  *>  v'{a>  ■ 


as  e  >  0  and 


,'<„,,e(x.t)>  =  -  ]'  *  ll  . 


where  n  is  outward  normal  to  ft  and 

<3 


f  -  y  in  ft 


in  ft-ft 

g 


Then  we  have 


l^x9x  ^yMy“g(x,t) 


^ +  4 

-1  ♦  qE(x,g£(x,t) ) 


-1  +  7  /q  xe(*€  -  y>dy 


Substituting  this  approximation  into  (6)  we  have  (10).  Mow,  the  second  term 
in  the  right  hand  side  of  (10)  is  called  the  integrated  penalty. 

3.2.  Discretization  of  the  penalized  problem 

The  penalized  problem  (8)  with  the  boundary  conditions  (3)  -  (5)  and  (9) 
is  discretized  by  the  finite  difference  method.  Also  the  free  boundary 
equation  (10)  is  solved  by  Euler's  method.  The  intervals  0  <  x  <  a  and 
0  <  y  <  b  are  divided  into  N  and  M  equal  subintervals  of  width  h.  The 
mesh  size  of  time  is  denoted  by  At.  We  use  the  following  notations  in  the 
discrete  systems 


-5 


-  ih, 
Yj  *  jh# 

\mkht 
*i,j,*  = 
9i,* 


°  <  i  <  N 
0  <  j  <  M 

m  J  m 

,  0  <  k 

^i'V 

xe(xi»yj#tk>  . 


We  define  the  discrete  characteristic  function  by 


i  -  p 


(ii) 


i.k 


.  h 

1  4e  pi,k 


j  >  tgi#k/hl 


j  -  [gi#k/h] 


j  <  Igi#k/h] 


where  [  ]  denotes  the  Gauss  symbol  and 

Pl,k  *  '  ’’i,!,"  • 

If  we  apply  five  points  formula  for  the  Laplacian,  we  have  the  following 
equations  the  potential  function  ♦  satisfies  for  any  k. 


<4  +  e  Xi, j.k^i, j,k  "  *i-1,j,k  "  *i+1,j,k 


(12) 


"  *i,j-1,k  “  *i,j+1,k  "  e  yj  Xi,j,k 


(0  <  i  <  W,  0  <  j  <  M)  . 


This  system  of  linear  equations  is  solved  by  the  incomplete  Cholesky 
decomposition  combined  with  conjugate  gradient  method  [13]. 

The  free  boundary  g^k  is  obtained  by 

"3)  ’i,*..  ■  ’x.k * “ ,(n.k>  <»<i<».  »<k>  . 


(14) 


xl,j,k'*itj,k  ‘  'j1 


It  should  b«  noted  that 


and 


are  determined  by  g^ 


3.3.  How  to  choose  the  penalty  parameter  e 


He  assume 

e  is  expressed  by 

(15) 

£  “ 

h° 

and  try  to  find 

an  optimal  value  of 

0 

right  sides  of  equations 

(6)  and  (14). 

test  problem: 

Au  m  0 

in 

(16) 

< 

u  *  1-y 

on 

u  =*  1-x 

on 

u  =  0 

V 

on 

This  problem  has 

an  exact  solution: 

u  “ 

1  ■ 

( a  >  0) 

to  minimize  the  difference  of  the 
For  this  purpose,  we  consider  a  simple 

ft 

g 

x  =  0 
y  =  0 
y  *  1-x  . 

x  -  y 


in  ft  (see  Figure  3). 


Figure  3 
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Therefore, 

(17) 


-  AT 


2  2u. 

9x  3n|y-g(x) 


where  g(x)  ■  1-x. 

e 

Here  we  construct  the  discretized  problem  (P  ) 

n 

equation  (PG ) : 

(18)  Au£  -  —  xu£  “  0  in  fl 

where 

( 1  in  a-a 

X-{ 

l  0  in  n 
v  g 


of  the  penalized 


We  investigate  the  difference  between  (17)  and  the  integrated  penalty: 


c  h  r 

qi  "  e  j£0  xi,jui,j 

by  varying  the  value  of  o  in  (15).  We  find  that  the  optimal  value  of  o 
3  ~  4  for  1/8  <  h  <  1/16. 


3.4.  Stability  condition 

Here  we  study  the  stability  condition  of  (13).  It  is  well  known  (14] 
that  the  stability  condition  is 

‘,9>  HUI  - 2  ■ 

where 

F(g)  “  +  ^  /q  xe(g)  (♦e(g)  -  y)dy  . 

If  we  use  the  property: 

I^U^lx^))!  <  C0/e  , 

where  CQ  is  a  constant  Independent  of  c  [15],  then  we  have 

cn 

I  ft  g)  -  F(g)l  <  —  Ig-gl  . 

/e 

If  we  substitut  Cfl//e  to  |dF/9g|  in  (19),  then  we  have 


is 


0  <  At  —  <2  . 

/  €  ” 


Therefore  we  may  choose 
(20) 


At  =  C^/e  *  C^h 


a/2 


4.  NUMERICAL  EXAMPLE 

In  this  section  we  present  results  for  the  problem  (1)  -  (6),  obtained  by 
the  method  of  Integrated  penalty.  Data  of  the  problem  are  as  follows: 

a  =  1 

b  -  0.3125 
A  -  0.625  . 

The  tapping  hole  is  located  at  (0,  1/16).  The  initial  surface  is  given  by 

g(x,0)  “  0.25 

The  parameters  used  for  the  numerical  calculations  are 

h  -  1/16 
e  “  h3  *  1/4096 
At  =  /c  “  1/64  . 

In  Figure  4  we  show  the  results. 


Figure  4 
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